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We construct a set of states that implement the non-Abelian Gauss’s law for QCD. We also construct a set of gauge-invariant 
operator-valued quark and gluon fields by establishing an explicit unitary equivalence between the Gauss’s law operator and the 
‘pure glue’ part of the Gauss’s law operator. This unitary equivalence enables us to use the ‘pure glue’ Gauss’s law operator to 
represent the entire Gauss’s law operator in a new representation. Since the quark field commutes with the ‘pure glue’ Gauss’s law 
operator, it is a gauge-invariant field in this new representation. We use the unitary equivalence of the new and the conventional 
representations to construct gauge-invariant quark and gluon fields in both representations, and to transform the QCD Hamiltonian in 
the temporal gauge so that it is expressed entirely in terms of gauge-invariant quantities. In that form, all interactions between quark 
fields mediated by ‘pure gauge’ components of the gluon field have been transformed away, and replaced by a nonlocal interaction 
between gauge-invariant color-charge densities. This feature — that, in gauge-invariant formulations, interactions mediated by pure 
gauge components of gauge fields have been replaced by nonlocal interactions — is shared by many gauge theories. In QED, the 
resulting nonlocal interaction is the Coulomb interaction, which is the Abelian analog of the QCD interaction we have identified 
and are describing in this work. The leading term, in a multipole expansion, of this nonlocal QCD interaction vanishes for quarks 
in color-singlet configurations, suggesting a dynamical origin for color confinement; higher order terms of this multipole expansion 
suggest a QCD mechanism for color transparency. We also show how, in an SU(2) model, this nonlocal interaction can be evaluated 
nonperturbatively. 


1 Introduction 

The program of the work we are discussing is de¬ 
signed to take advantage of a common thread that ap¬ 
pears to underlie all gauge theories — non-Abelian as 
well as Abelian: When only gauge-invariant fields are 
used in constructing the Hamiltonian for a gauge the¬ 
ory, interactions between charged fields (electrical, color, 
etc.) and pure-gauge components of gauge fields cannot 
arise, and a nonlocal interaction between charge den¬ 
sities — the Coulomb interaction in QED — appears 
in their stead. Moreover, in gauge-invariant QED, the 
only interaction besides the Coulomb interaction is be¬ 
tween the transverse (gauge-invariant) part of the gauge 
field and the current density; this interaction takes the 
form — fji(r)ATi(r)dr, where An is the transverse 
(gauge-invariant) part of the gauge field, and j('') = 
e^t(r) dip(r). Since the current density has a, vjc de¬ 
pendence in the nonrelativistic limit, the Coulomb in¬ 
teraction is, by far, the most important electrodynamic 
force in the low-energy regime. 

We will show that a very similar effect occurs in 
QCD. The main difference between the two cases is 
the fact that constructing gauge-invariant gauge fields 
is much more difficult in non-Abelian theories than in 
Abelian ones. Also, the nonlocal interaction in QCD 
is more complicated, and involves not only quark color- 
charge densities, but also gauge-invariant gluon fields. 
Significantly, in QCD, just as in QED, the only inter¬ 
action other than the nonlocal one involves the gauge- 
invariant gauge (gluon) field and a quark color-current 
density that is the non-Abelian analog of the current den¬ 
sity in QED, and that can be expected to also have a u/c 


dependence in the nonrelativistic limit. 

In this paper we will show how to implement the non- 
Abelian Gauss’s law, how to construct gauge-invariant 
quark and gluon states, and how to express the QCD 
Hamiltonian in terms of these gauge-invariant fields. We 
will also discuss the physical implications of the formal 
results we obtain. 


2 Implementing Gauss’s law for QCD 


The Gauss’s law operator in QCD is given by 
Din“(r)=‘pure-glue’ part 

g%v) =ftn“(r) +gr^M?(r)n=(r)+Jo“(r), (1) 

^ ^ 

JSC) 

where jo(r) = g'if\v)^tf:{v) is the quark color-charge 
density, and Jq (r) is the gluon color-charge density. To 
implement the non-Abelian Gauss’s law, we must con¬ 
struct states that are annihilated by the Gauss’s law op¬ 
erator, i.e. we must solve 

g“(r)|4') = 0. (2) 


In QED, it is easy to implement Gauss’s law, {9ini(r) -|- 
Jo(i’)} lb = 0, because ^^^(r) -I-jo(i’) and diU.^{v) 
are unitarily. equivalent, and ^^^(r) |^) = 0 is sim¬ 
ple to solveEl In QCD, this method for implementing 
Gauss’s law is not available. 0“(r) and 9in“(r) can¬ 
not be unitarily equivalent, because these two operator¬ 
valued quantities obey inequivalent commutator alge¬ 


bras; whereas [9in“(r),9in('(r')] = 0, 




igfo-bcgc^j.)g^j, — rO 
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We will initially implement the ‘pure glue’ form of 
Gauss’s law for QCD, 


An“(r)|'P) = 0, (3) 


by constructing a state |4') = for which 

L>in“(r)|\E') = {ain“(r) + Jq (r) j'l'|0) = 0; |0) repre¬ 
sents a state that ^ annihilated by 9^11“ (r) — the so- 
called ‘Fermi’ stateB We then seek to construct an oper¬ 
ator for which 


[a,n“(r), = (4) 


where BQ{r) is an operator that has 9^11“(r) on its ex¬ 
treme right. Eq.(^) is essentially an operator differential 
equation, in which the commutator, [9in“(r), d^], is a 
derivative. 


We have found a solution of this equation! 
is expressed as 


in which 


The condition on the resolvent gauge field that is 
equivalent to implementing the ‘pure glue’ Gauss’s law, 
is 


igr^<^A^{v) I dr'[Utir),^ir')]v;^{r') + 

1 1 dr'[d.Ut{r),^ir')]V;^{r') + A>^ir)Vf{r) 

= E f ) - 

r]=l 

2^2^l 9 r/!(t-l)!(i-|-l) ^ 

r/—0 t—1 

/(':?VS‘'Kf„(r)Mf„(r)a.r.Vr) 

r)—0 i—1 


4' = ||exp(y4) II 

(5) 

where B{g) is the Bernoulli number, and rep¬ 

resents any arbitrary vector field in the adjoint represen- 

with 

r 


tation of SU(3). 

A = i dr A] (r) B/ (r), 

(6) 

We have solved Eq.(||);i The solution is given by 

and with A] (r) represented as the series 


r OO n 

J dr A] (r)E/(r) = E ^ { ^G)j(’') + 

37(r) = - 

n—1 

(7) 

/,"fA<f„(r)E«„(r)}y7r), (10) 


The are nonlinear functionals of transverse and 

longitudinal parts of gauge fields, but are independent 
of the canonical momentum Aj (r). The ordered prod¬ 
uct II exp(yl) II is defined so that, in the order term, 
II (yl)” II , all functionals of the gauge field ^4“ are to the 
left of all functionals of the canonical momenta IlJ . We 
refer to Af (r) as the resolvent gauge field. 

The requirement that |'I') implement Gauss’s law can 
be translated into a condition on the resolvent gauge 
field. Before we formulate this condition, we first define 
the following quantities: 

^“(r) = |^kl“(r) and 7^“)(r) = X^m]^ 

which are functionals of gauge fields; 3^“(r) = ^yl“(r) 
and = Om^i which have structures 

similar to X°^ (r) and 7^“^^ (r) respectively, but which are 
functionals of the resolvent gauge fields. We also need to 
define: 


ja/37 _ ja[l]/3b[l] j:b[l]a[2]b[2] y6[2]a[3]h[3] ^ . 

^ jb[i)-2]a[77-l]b[77-l] jb[ri-l]a[ri]j . 

these are chains of structure constants whose ‘links’ are 
summed over repeated indices. 


where 


= (-1)" 2fr,)*W (11) 

+ (^2;f(r)] , (12) 

1 

~ (n+1) ~S^ ) ('') ; (13) 


af and xf designate the transverse and longitudinal parts 
of the gauge field Af{r) respectively. Eq. ( 0 ) calls for 
summations over the multiplicity index g, which labels 
the multiplicity of 3^“(r) factors in and the mul¬ 

tiplicity of ff“(r) factors in 7^“^^(r). This multiplicity of 
factors makes Eq.(lO) a nonlinear integral equation that 
specifies the resolvent gauge field Aj (r) recursively. The 
possibility that this nonlinear equation has multiple so¬ 
lutions, and that these multiple solutions correspond to 
different physical regimes, is interesting, but as yet un¬ 
explored. 


In our earlier work,l3 we proved that Eq. ( 0 ) is a so- 
lution of Eq.(^) — a proof that we have referred to as the 
‘fundamental theorem’ — and that the resolvent gauge 
field, therefore, is the operator-valued field required for 
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the implementation of the non-Abelian Gauss’s law. We 
will see, in the remainder of this paper, that the resol¬ 
vent gauge field plays a pivotal role in the construction 
of gauge-invariant fields and in the transformation of 
the QCD Hamiltonian to a functional of gauge-invariant 
fields. 

3 Gauge-invariant quark and gluon fields 

The observation that underlies the construction of gauge- 
invariant quark and gluon fields is that the Gauss’s law 
operator, ^“(r), and the ‘pure glue’ Gauss’ law opera¬ 
tor Q'Uv) = are unitarily equivalent. In earlier 

work,l3 we have shown that 

Q-[v)=UcQ%y)Uc\ (14) 

where Ue = exp(Co) exp(C) with Cq and C given by 

Co = i j dr X°‘(r) Jq {r) , and (15) 

C = zy dr>’“(r)jp (r) . (16) 

We are therefore free to interpret the ‘pure glue’ Gauss’s 
law operator t/“(r) as the complete Gauss’s law operator, 
(y“(r), in a different, unitarily equivalent representation. 
We will refer to the conventional representation, in which 
17“ (r) is the complete Gauss’s law operator and 17“ (r) the 
‘pure glue’ Gauss’s law operator, as the C representation; 
and the new representation, in which 17“ (r) is the unitar¬ 
ily transformed complete Gauss’s law operator, as the J\f 
representation. It is manifest that the spinor (quark) 
field '(/'(r) commutes with 17“(r). Since the Gauss’s law 
operator is the generator of gauge transformations, 
is a gauge-invariant spinor (quark) field in the J\f rep¬ 
resentation. To find the corresponding gauge-invariant 
quark field in the C representation, we apply the trans¬ 
formation that appears in Eq.(p^, and obtain 

V’Gi(r) =^cV’(r)^C ^ = Vc(r)V’(r) (17) 

where 

Vc(r) = exp(-f53^(r)^) exp(-igA’“(r)^) , (18) 

and the represent the Gell-Mann SU(3) matrices. 
Vc(r) could be written as exp{—igZ°‘{r)^), where the 
Baker-Hausdorff-Campbell theorem can be used to ex¬ 
press Z°‘ as a functional of Tf" and y°‘. Vc(r) takes 
the form of a unitary operator that carries out a gauge 
transformation on the spinor field ip(j), but, in this case, 
Z°‘ is not a c-number function in the adjoint repre¬ 
sentation of SU(3), but is a complicated functional of 
the gauge field. Under a gauge transformation, pre¬ 
cisely compensating transformations are made on ijj(r) 


and Vc(r), so that ^/’Gi(r) is strictly gauge-invariant. 
When they appear In the Af representation, the color- 
charge and color-current densities, jg (r) = gtp'^ {r)^tp{r) 
and j“(r) = g^p^r)ai^^|;{r) respectively, therefore are 
gauge-invariant, although, in the C representation, both 
of these quantities transform gauge-covariantly, as vec¬ 
tors in the adjoint representation of SU(3). In the JV 
representation, the quark field '0(r) implicitly includes 
enough of the gluon field to achieve this gauge invari¬ 
ance. 

We have used Eq. ( 0 ) to find the gauge-invariant 
gluon field 

[4,.(r) ^] = Ifc(r) [Atir)^]V,-\r)+^Vc{r)d.Vc-\r) , 

n (19) 

or, equivalently, □ 

A^„(r)=A^,(r) + [5.,-^]:^(r). (20) 

We observe that, as in QED, the gauge-invariant gauge 
field Aq ■ (r) is transverse. But it is not merely the trans¬ 
verse part of the gauge field. In contrast to the gauge- 
invariant gauge field in QED, AQj(r) also involves the 
transverse part of the resolvent gauge field. 

We have expanded iIjq\ (r) and Aq - (r), and have veri¬ 
fied that our gauge-invariant quark and gluon fields agree 
with the perturbative calculations of Lavelle, McMullan, 
et. al. to the highest orderJx) which their perturbative 
calculations were available. Bi3 

4 A gauge-invariant QCD Hamiltonian 

We have been able to express the QCD Hamiltonian 

Hqcv = / * { in(‘(j')n“(r) -H iE)“ (r)E)“(r) + 

[(3m - iai {d^ - igA^{r)^)] tp{r)} , (21) 

entir^ in terms of gauge-invariant quark and gluon 
fieldsp by systematically transforming it, term by term, 
from the C to the Af representation. Under this transfor¬ 
mation, tp{r)^V^^{r)tp{r), ijj^{r)^tp^r)Vc{r), but the 
gauge field remains untransformed. Shifting Vc(r) to 
the right and V^“^(r) to the left until they encounter 
the A“(r)^, turns the second line of Eq.(|^) into 
f drip^r) {(3m — iaidi) t(}{r) Hj-a, with 

f ' - ^—h -' 

Hj-A = - J dr gV’^(r)a,4V’(r) ^L(i’) • (22) 

F“ (r)F“ (r) remains untransformed in this process, but 
we found that this similarity transformation, when ap¬ 
plied to n“(r)n“(r), generates nonlocal interactions bfe 
tween color-charge densities which we will discuss belowB 
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The QCD Hamiltonian that results from the trans¬ 
formation to the M representation is 

Hqcv =Jdr{ ^n“(r)n“(r) + (r)F“ (r) + 

i^\r)[/3m - iaidi]'ip{r)} + H' , (23) 


gauge field Hq^ (x). Hg-q — the other nonlocal inter¬ 
action in Hlh — couples quark to gluon color charge 
density. In Hg_Q the quark color-charge density is cou¬ 
pled, through the same Green’s function to a 

gauge-invariant expr^^ion describing ‘glue’-color, which 
we have found to beB’EI 


where H' describes interactions involving the gauge- 
invariant quark field. The parts of H' relevant to the 
dynamics of quarks and gluons can be expressed as 


H' = Hg_A + Hlr , 

(24) 

where Hlr is the nonlocal interaction 


Hlr = Hg-q + Hq-q . 

(25) 


We will initially focus our attention on Hq-q in this 
report, since it illustrates the most important features of 
ou£| work. A useful formulation of Hq-q can be given 
asB 

Hq-Q = ^J drdxj(;(r)JC'^“(r,x)j(((x). (26) 

where the Green’s function !F^°‘ (r, x) is represented as 


jr'’“(r,x) 


Sab 

47r|r — x| 


K“(r)=5r"^Tr ^-^(r)VHc(r)^ Ag,,(r)n?(r). 


The chains that constitute x) — coupled links, 

in which the order chain is a product of n links 
— suggest features closely associated with QCD: flux 
tubes, ‘string’-like structures tying colored objects to 
each other, etc. But these rudimentary analogies only 
serve to direct our attention to potentially important fea¬ 
tures of this interaction — they do not, by themselves, 
demonstrate a physical effect. 

With regard to Hj-A: we observe that the color- 
current density j(‘(r) has a configuration-space structure 
that is similar to that of the electric current density in 
QED. We can therefore reasonably expect that it, too, 
will manifest a u/c dependence for nonrelativistic quarks 
interacting with the purely transverse Aq^ (x), and that 
Hlr will be of predominant importance in the low-energy 
QCD regime. 


2g 


dy 


47r|r — y 




1 


47r|y-x| 




4,‘.’(yi)a 


dy2 


47r|r - yi 


47r|yi - y2| 


4r’(y2)9, 


dys 


dyn 


^Q\e. 


47r|y2 - ya 
(y»)9( 


47r|y(„_i) - yr, 


1 


47r|y„ -x| 


(27) 


We make the following observations about x): The 

initial term resembles the Coulomb interaction. The in¬ 
finite series of further terms consists of chains through 
which the interaction is transmitted from one color- 
charge density to the other. Each chain contains a suc¬ 
cession of ‘links’, which have the characteristic form 


5 Quark confinement and color transparency 


The use of Eqs.(p^ and (|^ to explicitly evaluate an 
effective ‘potential’ between quark color-charge densities 
is precluded, at the present time, by the fact that we have 
not yet found an expression for the gauge-invariant gauge 
field Aq^(x). However, we can draw some important 
conclusions from the general form of Eqs.(26) and (HI)- 


If we assume that the Green’s function .F^“(r, x) 
varies smoothly as a function of x, when evaluated in 
an appropriately chosen state, and that quarks, or con¬ 
figurations of quarks, are localized in wave packets whose 
size is small enough so that .7^^“(r, x) varies only moder¬ 
ately over the space occupied by these wave packets, then 
we are entitled to make a ‘color-multipole’ expansion of 
the expression /lF^“(r, x)jQ (x)(ix — which appears in 
Eq.(|^) — about the point x = xq, in the form 



{^^^■^“(r, xo) -h Xi5i.F'’“(r, xo)-h 


link = grw^^(=)A^i,(x)a 


47r|x — y| 


(28) 





(29) 


The ‘links’ are coupled through summations over the 
S(„) indices and integrations over the spatial variables. 
As noted before, all the quantities in Hq-q are gauge- 
invariant — the color-charge density Jq (r) as well as the 


where Xi = {x — a:o)i and di = d/dxi. When we perform 
the integration in Eq. ( |2^ , the first term contributes 
_^f)o(r, xp) , where Q'^ = f dxjg(x) (the integrated 
“color charge”). Since the color charge is the generator 
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of rotations in SU(3) space, it will annihilate any multi¬ 
quark state vector in a singlet color configuration. Multi¬ 
quark packets in a singlet color configuration therefore 
are immune to the leading term of the nonlocal Hq-q. 
Color-singlet configurations of quarks are only subject 
to the color multipole terms, which act as color analogs 
to the Van der Waals interaction. The scenario that this 
model suggests is that the leading term in Hq-q, namely 
Xo)Q“ for a quark color charge Q“ at Tq and 
another quark color charge at Xq, may be responsi¬ 
ble for the confinement of quarks and packets of quarks 
that are not in color-singlet configurations. Moreover, 
assuming that Xq) varies only gradually within 

a volume occupied by quark packets, the effect of the 
higher order color multipole forces on a packet of quarks 
in a color-singlet configuration becomes more significant 
as the packet increases in size. As small quark pack¬ 
ets move through gluonic matter, they will experience 
only insignificant effects from the multipole contributions 
to Hq_q, since, as can be seen from Eq. (29), the fac¬ 
tors Xj , XjX. 


0 ■> 


keep the higher order 


multipole terms from making significant contributions to 
J’(ix.A^“(r, xjjQ (x) when they are integrated over small 
packets of quarks. As the size of the quark packets in¬ 
creases, the regions over which the multipoles are inte¬ 
grated also increases, and the effect of the multipole in¬ 
teractions on the color-singlet packets can become larger. 
This dependence on packet size of the final-state inter¬ 
actions experienced by color-singlet states — i. e. the in¬ 
creasing importance of final-state interactions as color- 
singlet packets grow in size — is in qualitative agreement 
with the characterizations of color transparency given by 
Miller and by Jain, Pire and Ralston. Ej 


In making this argument, the following features of 
this gauge-invariant formulation must be borne in mind: 
The state vectors that can appropriately be chosen for 
such a representation must implement the non-Abelian 
Gauss’s law applicable to QCD — i. e. the states would 
have to be annihilated by the Gauss’s law operator, which 
is Q in the N representation. That still allows us to apply 
quark creation operators to such states to form multi¬ 
quark Fock states with orbitals whose configuration- 
space dependence is arbitrary as far as the gauge problem 
is concerned. These states are far from being Fock states 
in the gluon sector — they have to be complex enough 
to implement Gauss’s law. But that fact does not inhibit 
us in using Fock states for multi-quark configurations, 
because, in the M representation, the quark fields con¬ 
tain enough gluon contributions to be gauge-invariant by 
themselves. The multi-quark states therefore have a sta¬ 
tus similar to multi-electron states in a gauge-invariant 
formulation of QED {e.g. in the Coulomb gauge), which 
are surrounded by a Coulomb field and therefore are 


gauge-invariant. 

In order to test these ideas quantitatively, we have 
made use of Yang-Mills theory — the SU(2) version of 
this model — for which the structure constants are 
We model Aq-( r), which, in this case, are trans¬ 
verse fields in the adjoint representation of SU(2), as the 
manifestly transverse “hedgehog” configuration 

(^GuW) = (30) 

Although there is no reason to believe that the ansatz 
given in Eq. (|^ ) follows from the dynamical equations 
that determine AQj(r), it is a convenient choice for ex¬ 
amining to what extent the structure of x) — as 

distinct from the precise form of AQj(r) — enables us 
to nonperturbatively evaluate the infinite series given in 
Eq. As we were able to showp the perturbative 

evaluation of Hq-q for this SU(2) model could be by¬ 
passed, and the entire series in X^°'(r,x) obtained from 
the solution of a sixth-order differential equation. 

Acknowledgments 

This research was supported by the Department of En¬ 
ergy under Grant No. DE-FG02-92ER40716.00. 


References 


1. P. A. M. Dirac, Can. J. Phys. 33, 650 (1955); 

K. Haller and E. Lim-Lombridas, Found, of Phys. 
24, 217 (1994), and further references cited therein. 

2. for example, R. Jackiw, Rev. Mod. Phys. 52, 661 
(1980); K. Haller, Phys. Rev. D36, 1839 (1987). 

3. E. Fermi, Rev. Mod. Phys. 4, 661 (1932). 

4. M. Belloni, L. Chen and K. Haller, Phys. Lett. 
B373, 185 (1996). 

5. L. Chen, M. Belloni and K. Haller, Phys. Rev. D55, 
2347 (1997). 

6. M. Lavelle and D. McMullan, Phys. Lett. B329 
(1994) 68. 

7. E. Bagan, M. Lavelle, B. Fiol, N. Roy and D. Mc¬ 
Mullan, Constituent Quarks from QCD: Perturba¬ 
tion Theory and the Infrared. hep-ph/9609330| . 

8. M. Belloni, L. Chen and K. Haller, Phys. Lett. 
B403, 316 (1997). 

9. L. Chen and K. Haller, Quark confinement and 
color transpa rency in a gaug e-invariant formula¬ 
tion of QCD, |rep-th/980325C . 

10. G. A. Miller, Color Transparency — Color Coher¬ 
ent Effects in Nuclear Physics. DOE/ER/41014- 
26-N97; nucl-th/970704C ; P. Jain, B. Pire, and J. 
P. Ralston, Physics Reports 271, 67 (1996). 


5 










